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ABSTRACT 



The radioactive isotope ^^Na is one of the primary 
tracer isotopes that may potentially be detectable 
by gamma-ray spectroscopy. This isotope is pre- 
dicted to be produced in ONe nova explosions, yet 
no detection of its 1.275 MeV gamma-ray signature 
has been reported so far. In particular, 6 years of 
COMPTEL observations of the galactic bulge, a re- 
gion that should be particularly rich in nova explo- 
sions, have only led to an upper-limit. In this pa- 
per we will present the implications of the COMP- 
TEL upper flux limit on the permissible ^^Na yield 
for ONe novae. We base our analysis on extensive 
Monte-Carlo simulations of galactic 1.275 MeV emis- 
sion that were constrained by the most recent results 
of galactic nova rates and distributions. We demon- 
strate that the non-detection of the 1.275 MeV line 
from the galactic bulge implies a solid upper ^^Na 
yield limit in agreement with current theoretical nu- 
cleosynthesis calculations. 
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1. INTRODUCTION 



Classical ONe nova outbursts are expected to synthe- 
size ^^Na and provide a 7-ray signal at 1.275 MeV 
(Clayton & Hoyle (1974), Gomez-Gomar et al. 
(1998)). Observations with several instruments have 
reported only upper-limits on the 1.275 MeV flux 
from the Galaxy (HEA03, SMM, OSSE, COMP- 
TEL) or from individual novae (COMPTEL) and 
consequently allow to derive an upper-limit of the 
^^Na mass ejected per ONe novae. Observations with 
the HEAOS 7-ray spectrometer provided an upper- 
limit of 4 10~^ photons cm~^ s~^ on the Galactic 
accumulated flux leading to an upper-limit of the 
ejected ^^Na mass of 3 10^*" Mq (Higdon & Fowler, 
1987). Leising et al. (1998) found a 99% confi- 
dence limit of 1.2 10~^ photons cm~^ s~^ on a steady 
1.275 MeV flux from the Galactic center (GC) direc- 



tion, lyudin et al. (1995), using COMPTEL obser- 
vations of single novae, estimated a 2a upper-limit of 
3 10~^ photons cm~^ for any neon-type novae in 
the galactic disk, which has been translated into an 
upper-limit of the ejected ^^Na mass of 3.7 10^^ Mq. 
Using observation of the GC with OSSE, Harris 
(1997) estimated an upper-limit of the ^^Na mass 
of 2 IQ-'^ Mq. Recently, lyudin et al. (1999) de- 
rived a 2a upper-limit of the ejected ^^Na mass of 
3.6 10-9 Mo and 2.1 10"^ M© from the COMPTEL 
observations of the bulge and Nova Cygni 1992, re- 
spectively. However, most of the ^^Na mass upper- 
limits derived from the observation of the Galactic 
diffuse emission have been calculated using overesti- 
mated ONe nova frequencies (up to 40 ONe novae 
per year only in the bulge). Current estimations of 
the ONe nova frequency range from 3 yr""'^ to 15 yr"^ 
in the whole Galaxy (see details of these figures in 
section 2). Moreover, the flux at 1.275 MeV in the 
GC region depends not only on the amount of ^^Na 
ejected per outburst and the rate of Galactic ONe 
novae but also on their spatial distribution. 

In the presented work, we estimate the upper-limit 
of the ejected ^^Na mass derived from COMPTEL 
observation of the cumulative 1.275 MeV emission in 
the GC region taking into account recent results of 
rates and spatial distributions of novae in our Galaxy. 
We modeled the Galactic emission at 1.275 MeV as 
a function of the ONe novae rate and the mean ^^Na 
yield per outburst. This has been done for currently 
used spatial distributions of novae in the Galaxy. 
The upper-limit of the ejected ^^Na mass is esti- 
mated as a function of the Galactic ONe nova fre- 
quency. The next section describes the method used 
to model the Galactic emission at 1.275 MeV. The 
estimation of the ^^Na mass upper-limit derived from 
the COMPTEL measurement is presented in section 
3. Discussion and conclusions follow (section 4). 



2. MODELING OF GALACTIC 1.275 MEV 
EMISSION 



We have used a method based on the work from 
Higdon & Fowler (1987) which estimated the first 
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Table 1. Characteristics of the models adopted for 
the simulations of the emission at 1.275 MeV from 
galactic ONe novae. R1/2 is the half radius of the 
spheroid, pt the radial scalelength of the disk, Zh the 
scaleheight of the disk and pb the proportion of novae 
that occur in the bulge. 



Models 


Ri/2(kpc) 


Ph{kpc) 


Zh{kpc) 


Pb 


HF87 


2.7 


3.5 


0.106 


0.348 


VdK90 


2.7 


5.0 


0.300 


0.105 


KDF91 


1.4 


3.0 


0.170 


0.179 


DJ94 


1.4 


5.0 


0.350 


0.111 



^^Na mass upper-limit using observation of diffuse 
1.275 MeV emission with HEA03. The distribution 
of the 1.275 MeV emission from galactic ONe novae 
is simulated with a Monte-Carlo method. The po- 
sition in galacto-centric coordinates and the age of 
ONe novae are chosen randomly according to the ap- 
propriate distributions and rates. For a given value 
of the ^^Na yield per outburst, the 1.275 MeV flux 
arriving on Earth from each nova is calculated with 
the distance and the age of the nova (see Jean et al. 
(2000) for a detailed description of the method). 

For the purpose of this work, a Galaxy model with 
a disk and a spheroid is convenient to simulate the 
distribution of Galactic ONe nova events. Several 
laws for the spatial distribution of classical novae 
in the disk and in the 'spheroid' (representative of 
the bulge) have been proposed. We have selected 
four models that differ signiflcantly from each other. 
The first and older of them is described in Higdon 
& Fowler, 1987 - hereafter HF87. The disk and 
spheroid models are derived from the starlight sur- 
face brightness distribution. By scaling the nova rate 
in the bulge of M31 to the bulge of our Galaxy, HF87 
estimate the proportion of novae in the spheroid to 
0.348. The second model has been used by Hatano et 
al. (1997) to estimate the occurrence rate of Galac- 
tic classical novae and the fraction of novae in the 
bulge. It is based on a model of the distribution of 
type la supernovae by Dawson & Johnson (1994) - 
hereafter DJ94. The proportion of novae that occur 
in the bulge is set to 0.111 on the basis of an es- 
timate of the bulge to total galaxy mass ratio. The 
third model is derived by Kent, Dame & Fazio (1991) 
- hereafter KDF91 - from the Galactic survey of the 
Spacelab InfraRed Telescope (IRT) that provides a 
reliable tracer of the distribution of G and K giant 
stars. Using the total infrared luminosity of the bulge 
and the disk, the derived proportion of novae occur- 
ring in the bulge is 0.179. The last model is taken 
from Van der Kruit (1990) - hereafter VdK90. It has 
been used by Shafter(1997) to estimate the galactic 
nova rate. This author assumes that the nova distri- 
bution follows the brightness profile of our Galaxy. 
Under this assumption, the proportion of bulge no- 
vae is 0.105. Table |l| summarizes the characteristics 
of the adopted models. 

The Galactic nova rate is poorly known, indepen- 
dently of the underlying white dwarf composition, 



because the interstellar extinction prevents us from 
directly observing more than a small fraction of the 
novae that explode each year. Several methods have 
been used to estimate the nova occurrence rate. Esti- 
mations based on extrapolations of Galactic nova ob- 
servations suggest a rate in the range 50-100 yr~^(see 
Shaffer (1997) and references therein). Observations 
of novae in other galaxies have revealed a correlation 
between the nova rate and the infrared luminosity 
of the parent galaxy. However, the occurrence rates 
of Galactic novae computed by scaling the nova rate 
measured in external galaxies (« 20 yr^^, see Daw- 
son & Johnson (1994)) are lower than 50 yr~^. Re- 
cently, Shafter (1997) reconciled this difference by 
recomputing the nova rate with the Galactic nova 
data. He extrapolated the global nova rate from the 
observed one, accounting for surface brightnesses of 
the bulge and the disk components and correction 
factors taking care of any observational incomplete- 
ness. With this method, Shafter (1997) estimated 
the nova ratetobe35±llyr-^ Hatano et al. (1997) 
found a similar value (41±20 yr~^) using a Monte- 
Carlo technique with a simple model for the distri- 
bution of dust and classical novae in the Galaxy. 

The Galactic ONe nova rate is obtained by mul- 
tiplying the total nova rate by the fraction of no- 
vae that results from thermonuclear runaways in 
accreted hydrogen-rich envelopes on an ONe white 
dwarf. Several authors deduced from observations of 
abundances in nova ejecta a proportion of ONe no- 
vae from 20 to 57 percent of observed nova outbursts. 
Livio & Truran (1994) reestimated the frequency of 
occurrence of ONe novae, in light of observations of 
abundances in nova ejecta. They concluded that, of 
the 18 classical novae for which detailed abundance 
analyses were available, only two or three had a large 
amount of neon and were ONe novae, whereas three 
other novae showed a modest enrichment in neon, 
casting doubt on the type of the underlying white 
dwarf. Under these considerations, they estimated a 
fraction of Galactic ONe novae between 11 and 33 
percent. 

Therefore, for the calculations presented in the next 
section, we have adopted frequencies of Galactic no- 
vae ranging from 24 yr~^ to 46 yr~^ and a proportion 
of ONe novae from 11 % to 33 %, in agreement with 
the recent estimates of Livio & Truran (1994) and 
Shafter (1997). These values correspond to a lower 
and upper- limit of the ONe nova rate of 3 yr~^ and 
15 yr~^, respectively. Figures |l| show examples of 
computed 1.275 MeV line fiux spatial distribution 
for the four chosen models and a galactic ONe nova 
rate of 10 yr^^. 



3. ESTIMATION OF THE ^^NA UPPER-LIMIT 

Several Monte-Carlo simulations have been done for 
a given galactic ONe-novae frequency-spatial distri- 
bution and a fixed value of the ^^Na yield per nova. 
For a large number N^ of Galaxy-tests, the prob- 
ability for a detection of the cumulative emission 
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Figure 1. Examples of modeled diffuse 1.275 MeV 
emission generated by Monte Carlo simulation for 
the four chosen spatial distributions. The adopted 
ONe nova rate is 10 yr^^ . The fluxes have been nor- 
malized with the '^^Na ejected mass. The intensity 
is plotted with 3° by 3° pixels in a map in galactic 
coordinates. 
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Upper limit of the ''Na mass ejected per outburst (M^ ) 

Figure 2. Relation between the "^"^Na mass upper-limit 
derived from the CO MP TEL 1.275 MeV flux upper- 
limit in the ICP radius around the GC, and the Calac- 
tic ONe nova rate, for the 4 adopted Galaxy models. 



has been estimated by calculating the number of 
Galaxy-test for which the simulated flux in a circle 
of 10° radius around the GC is above the upper-limit 
flux of COMPTEL (/2<T=3.0 IQ-^ photons cm^^ g-i 
- lyudin et al., (1999)) during the lifetime (6 years) 
of the instrument. The probability is obtained by di- 
viding Ni with the total number of Galaxy-test Nj-. 
Such fractions above f2a can be interpreted simply 
as the fraction of time COMPTEL would detect such 
flux. The ^^Na mass upper-limit is defined when 
the probability of detection of the 1.275 MeV flux 
by COMPTEL reaches 90%. It has been calculated 
for several galactic ONe nova rate and for the pro- 
posed spatial distribution models. The results are 
presented in the figure ||. 



4. DISCUSSION 



The ^^Na mass upper-limit derived from the 
1.275 MeV emission of the bulge depends strongly 
on the adopted ONe nova rate and spatial distribu- 
tion. Based on recent results on ONe nova rate, we 
show that this upper-limit can not be lower than 
2 10~^ M0 per nova (see figure ^). Moreover, tak- 
ing into account the uncertainty in the ONe nova 
frequency in our Galaxy, the most reliable upper- 
limit should be estimated using the lowest ONe nova 
rate value and the less favorable spatial-distribution. 
Consequently, with an ONe nova rate of 3 yr~^ and 
the VdK90 model we obtain a ^^Na mass upper- 
limit of 6 IQ-'^ Mq with the COMPTEL measure- 
ment. Indeed, if the yield of ^■^Na is larger than 
this upper-hmit, COMPTEL should have detected 
the 1.275 MeV line in the bulge especially if the ONe 
nova rate is larger and if the distribution is more fa- 
vorable for the detection of this 7-ray line. Applying 
the same analysis with the 1.275 MeV flux upper- 
limit derived by the observation of the central 60° 
of the Galaxy by the IIEA03 spectrometer (Higdon 
& Fowler (1987)) yields a ^^Na mass upper-limit of 
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Figure 3. Relation between the Na mass upper- 
limit, derived from the HE AOS Galactic diffuse 
1.275 MeV fl,ux upper-limit (in a longitudinal band 
ranging from 3 3 CP to 3 CP ), and the Galactic ONe 
nova rate, for the 4 adopted Galaxy models. The dot 
shows the estimation of Higdon & Fowler (1987). 

2.5 10-6 Mq instead of 3 10"^ Mq (see figure |). 

The ^^Na mass upper-limit derived by lyudin et al. 
(1999) is «170 times lower than the value obtained 
from our analysis of COMPTEL measurements. The 
difference can be explained by the different 0-Ne 
nova rate that has been assumed for the galactic 
bulge region. While lyudin et al. (1999) assumed 40 
0-Ne novae per year in the bulge region, we estimate 
this rate to only about 0.3 per year (see table 1 and 
section 2). lyudin et al. (1999) simply employed the 
'canonical frequency of 40 novae yr~^' in the Galaxy 
as the 0-Ne nova frequency of the galactic bulge re- 
gion. Ho"wever, only about 1/4 of all galactic novae 
are of 0-Ne type, and also, only a small fraction of 
them is located within the galactic bulge. Hence we 
believe that our more conservative estimate is more 
realistic, and the actual constraints imposed by bulge 
gamma-ray observations on nova ^^Na yields are less 
stringent than claimed by lyudin et al. (1999). 

Since the pioneering work of Clayton & Hoyle (1974), 
pointing out the potential role of novae as gamma-ray 
emitters, several theoretical nucleosynthesis studies 
have been performed to predict the amount of ^^Na 
ejected in nova outbursts. They have shown that 
only ONe novae are likely sites for the synthesis 
of ^^Na due to the higher peak temperatures at- 
tained during the explosion (with respect to CO 
novae) but also due to the presence of seed nuclei 
mainly ^"Ne, extremely relevant to ^^Na synthesis 
since CNO breakout scarcely takes place at typi- 
cal nova temperatures. Recent fully hydrodynami- 
cal calculations of nova outbursts with updated nu- 
clear reaction rates indicate that nova release about 
7 10-9 (for 1.15 Mq white dwarf) to 4.4 lO'^ Mq 
(1.35 Mq white dwarf) of ^^Na into the interstel- 
lar medium during the explosion (Jose & Hernanz 
(1998), Jose, Coc & Hernanz (1999)). Such values 
are below the limits posed by several attempts to 
find the 1.275 MeV gamma-ray signature from clas- 



sical novae. Although nuclear uncertainties associ- 
ated to key reactions relevant to ^^Na synthesis ex- 
ist (i.e., ^^Na(p,7) and ^^Na(p,7)) their incidence in 
the expected yields is small, namely a factor of 3. 
A recent reanalysis of the uncertainties associated 
to 2iNa(p,7) (Smirnova & Coc (2000), Coc et al. 
(2000)) suggest even a smaller impact on the ^^Na 
yields. 

The most stringent ^^Na mass upper-limit obtained 
so far is presented by lyudin et al. (1999) who de- 
rived a 2(7 upper-hmit of 2.1 10"^ M© from COMP- 
TEL observations of Nova Cygni 1992. This upper- 
limit still remains above the most optimistic the- 
oretical estimates. Observation with future 7-ray 
spectrometers should provide more stringent upper- 
limit. Jean et al. (2000) show that the likelihood 
to detect the Galactic diffuse 1.275 MeV emission 
with SPI (Spectrometer of the INTEGRAL mission) 
with only wlO days of observation is small. On the 
other hand, they estimate that 80 days of observation 
of the GC region could already provide constraints 
to the mean ^^Na yield in ONe novae as a function 
of their frequency-spatial distribution in the Galaxy. 
Using the most recent calculations of ^^Na yields in 
ONe novae and SPI instrumental background pre- 
dictions, Hernanz et al. (2000) estimate that the 
1.275 MeV line from an individual ONe nova could 
be detected by the spectrometer of INTEGRAL if 
its distance is less than 1 to 2 kpc, depending on the 
uncertainties of nuclear reaction rates. 

Gamma-ray observation of novae would provide in- 
formation not only on their eruption mechanisms and 
the nucleosynthesis processes involved in their explo- 
sion but also on their distribution and their rate in 
the Galaxy (e.g. proportion in the bulge, scaleheight 
in the Galactic disk) since the problem of the in- 
terstellar extinction does not appear at this energy 
range. 
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